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taking into account dead-cone effects and soft radiation characteristic of massive charges. 
This factor multiplies the massless jet distribution function Jn(Q 2 )- In the above equation 
the variable iV is rescaled by the mass correction parameter m 2 /Q 2 <C 1 with respect to the 
standard massless case, being m the quark mass and Q the hard scale. The functions A(as) 
and B(as) appear with a minus sign suppressing collinear effects at very large N > Q 2 /m 2 , 
as expected. In the same region, soft radiation not collinearly enhanced, characteristic of 
on-shell massive charges, makes its appearance with the function D(as)- Phenomenological 
applications, such as the resummation of b — > civ decay spectra or the inclusion of beauty 
mass effects in t — > bW decays, are briefly sketched. 



* This work is dedicated to the memory of Jiro Kodaira. 



1 Introduction 



dP~ a —, (2) 



The structure of hadronic final states at high energy is largely determined by the infrared (soft 
and collinear) divergencies occurring in real contributions to QCD cross sections. The singular- 
ities cancel with those coming from virtual corrections in sufficiently inclusive observables, but 
leave often large residual logarithmic effects [1]. Soft singularities mean in physical terms a high 
probability of soft gluon emission: 

dP~*§, (1) 

where a = as with E the gluon energy and Q the hard scale (Aq C d <C E <C Q, Aq C d being 
the hadronic scale). Because of such singularities, one has to think to a quark or a gluon as a 
"dressed" parton, which is always accompanied by a cloud of virtual soft gluons. In a scattering 
process, part of the cloud "detaches" and converts itself into real soft radiation, observed as 
low-energy hadrons. Collinear singularities in physical terms mean high probability of collinear 
configurations: 

dff 2 

e 2 

where is the emission angle. Because of them, a massless quark or a gluon evolves into an 
ensemble of collinear partons, which show up in the detectors as a jet of hadrons. In the case 
of a massless parton, there is an overlap of soft and collinear regions producing the well-known 
(leading) double-logarithmic effects: 

dE dO 2 

Furthermore, there is hard collinear emission related to contributions of the form (2). In the 
case of a massive parton — which we assume to be observed in a reference frame where it is 
fast moving — two main changes occur in the QCD radiation. The first is the well-known "dead 
cone" effect, according to which little radiation is emitted in the forward direction: 

dP ^ % 2 + m 2 /Q 2 ~ aH { 9 -Q*) IT. W 

where H is the step function. The second effect is related to the fact that a massive parton 
radiates soft quanta isotropically in its rest frame: 

dP m ~ a ^r- ( 5 ) 

These quanta, being radiated in any space direction, cannot be ascribed in a natural way to 
any jet in the event. They are related to the "classical" chromo-electric field of a static color 
charge. By means of a Lorentz transformation, we obtain a fast-moving parton, which is now 
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accompanied by a boosted Coulomb field, i.e. by soft radiation restricted to a forward cone 
of angular size #~l/7£~m/Q. As a consequence, when we consider the structure of a 
jet initiated by a parton with a small mass, we expect missing collinear radiation compared 
to the massless case and additional soft radiation related to the boosted Coulomb field. A 
non-vanishing quark mass, m 7^ 0, produces, as we are going to show in detail, a specific sub- 
structure, which is expected to be universal on the basis of physical intuition: small angle 
emissions only are involved, which can be ascribed to a specific jet in the event. In general, 
whenever logarithmically enhanced effects are encountered, one expects a factorized structure. 
As we are going to show in detail in this note, all these expectations turn out to be correct. 

In the massless case, it is convenient to introduce the "jet function" J(y;Q 2 ), which gives 
the probability that a massless parton produced in a hard process characterized by the scale Q 
fragments into a hadronic jet of mass squared [2-4] 

m\ = yQ 2 . (6) 

In first order, its expression reads: 



J(y; Q 2 ) = 5{y) - A, a + B x a + O (a 2 ) , 



(7) 



where A 1 and B 1 are coefficients whose explicit expressions will be given later and the plus 
distributions are defined as: 



P(y)+ = lim 



9{y - e)P(y) - S(y - e) / P(y')dy 



(8) 



Note that in the free limit, a — > 0, the jet has a trivial structure and the mass distribution 
reduces to a spike corresponding to the (zero) parton mass. If y <C 1, the coefficients of a in 
eq. (7) become large, making a truncated perturbative expansion unreliable and asking for a 
resummation to all orders. The latter is systematically performed by going to A— space, i.e. by 
considering 

Jn(Q 2 ) = f dyil-yY-'Jiy-Q 2 ). (9) 
Jo 

The resummed expression reads [2-4]: 

Jn(Q 2 ) = exp / V 1 ~ ^ ~ * ( r Vd $A[a{kl)} + B [a (Q 2 y)] \. (10) 

Jo y [ Jq 2 v 2 ft ± j 

The functions A(a) and B(a), generalizing the first-order expressions A\a and B\ol respectively, 
describe small- angle soft and hard gluon emission off a massless parton respectively (see later) . 
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In the massive case, one can define a generalized jet function J(y; Q 2 ,m 2 ), where y is now 
defined as (cfr. eq. (6)): 



m x — m 



y = (») 



To first order, the massive jet function reads: 

6{y-r) 



J(y; Q\m 2 ) = S(y) + a 



{-A 1 logy + Bt) 



+ (-Ax logr + 



y 

6{r-y) 



y 



+ O {a 2 ) , (12) 



where is the step function and D 1 is a coefficient which will be derived in the next section. We 
have defined the mass correction parameter 

r S g « 1. (13, 

Note that r is quadratic in m and that we always assume that the quark mass is much smaller 
than the hard scale, m < Q, is in order to have fast-moving charges and to preserve a jet 
structure 1 . 

According to eq. (12), we can identify two different kinematical regions, which can be defined 
with logarithmic accuracy: 

1. high jet mass: y ^> r or, equivalently, mx —m^>m. The quark mass m can be neglected 
and the collinear region produces logy terms; 

2. low jet mass: ?/ C r or, equivalently, mx —m<^m. The quark mass screens the collinear 
singularity and produces logr terms. 

Mass effects are more easily looked at by considering the partially integrated jet rate: 

R(y; Q\ m 2 ) = f J(y>; Q 2 , m 2 ) dy' . (14) 
Jo 

In the high-mass region: 

fi 1 

R(y) = 1 - J J(y'; Q 2 ,m 2 )dy' = 1 - -A.alog'y + B 1 a logy (y > r), (15) 



x We are not interested to the case in which the quark mass is close to the hard scale, m w 0(Q), as this case 
corresponds to a non-relativistic motion of the color charges which does not give rise to jets. 
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while in the low-mass region: 



R(y) = 1 — Ai a \ogy \ogr + ^ a \og 2 r + Dialogy + (Bi — Di)a\ogr (y < r). (16) 

For r <C 1 mass logarithms become large and need to be resummed to all orders. We will show 
that mass terms can be relegated into a factor which takes into account the effects discussed 
above. The jet function can be factorized in moment space as: 

Jn(Q 2 ; m 2 ) = J N (Q 2 ) 5 N (Q 2 ; m 2 ), (17) 

where the mass-correction factor 5n(Q 2 ', rn 2 ) reads: 

f 1 (\ _ n) r(JV-l) _ I f rm 2 y ,.2 

6 N (Q 2 ;m 2 ) = exp / dy { -^> - - [a (k*)] — B [a (m 2 y)} + 

JO V [ Jm 2 y 2 fc _L 

+ D[a(m 2 y 2 )] J- (18) 

The function D(a), generalizing the first-order term Dia, describes soft radiation not collinearly 
enhanced off a nearly on-shell massive quark (see later). We derive our main result, eq. (18), with 
some rigour in next-to-leading logarithmic (NLL) approximation and we conjecture its general 
validity beyond NLL. Eq. (18) has a simple physical interpretation: the parameter N — 1 is 
multiplied by r on the r.h.s., implying that mass effects become "visible" only for large 

O 2 

N > ^ > 1. (19) 

In this case, there is enough resolution to "see" the quark mass, which tends to suppress the 
collinear effects, related to the A and B terms. At the same time, soft radiation not collinearly 
enhanced, described by the function D and characteristic of massive partons, does appear. Let 
us also note that, since the jet mass is an infrared (i.e. soft and collinear) safe quantity, 5n — 1 
for r = 0. 

For clarity's sake, we derive the general formula (18) in sec. 2 for a specific physical process 
first, namely the radiative decay 

B -> X, + 7 (20) 

by keeping m s ^ 2 . In sec. 3 we generalize the results of the previous section to the case 
of a process involving one or more jets initiated by quarks with a small mass compared to the 
hard scale. Finally in sec. 4 we present the conclusions of our study and we discuss natural 
applications and developments. 



2 Lct us note that the inclusion of strange mass effects is a rather academical problem because even with a 
large value of the constituent mass m s — 0.5 GeV, the correction parameter is very small: r — (m s /mb) 2 w 10~ 2 . 
Furthermore, strange mass effects become visible for x 7 = 2E 1 /rriB = 1 — m 'xs/' m % ~ 0-99, where resonance 
effects due to K, K* peaks, etc. — intrinsically beyond perturbative QCD — are substantial. 
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2 Mass Distribution in B — > X s + 7 for m s 7^ 



In this section we compute the invariant hadron mass distribution in the decay (20), having the 
fixed-order expansion: 

ldF . 1 dT^ 1 dT^ 

where 1/T dT^/dy = 0(a n ) and T is the corrected inclusive width. The massless case has been 
considered for example in ref. [5]. 



2.1 Single-gluon distribution 

In order to include the quark-mass corrections to eq. (21), let us start from the eikonal current 
containing explicitly the mass terms. According to standard factorization formulae, the soft 
terms in the amplitude for 

b -> s + g + 7 (22) 

are factorized by the eikonal current 

■™ = *(%&- T -i&)- (23) 

where g is the coupling in the QCD lagrangian, is the soft gluon momentum and T b and T s 
are the color generators for the beauty and the strange quark respectively [6]. Since the photon 
- in general the probe — carries no color, we have color conservation along the quark line: 

T b = T s . (24) 

The square of the eikonal current reads: 3 



g^J»(k) J v {k)* = ArraC F 



2 (p b ■ p 8 ) ml m 2 s 



(p b ■ k) (p 3 ■ k) {p h . k) [p s ■ ky 



(25) 



, C F f 1 + r 1 r 



E 2 g [1-r t + r/(l-r) (1 - r) 2 [ t + r /(l - r)Y 

where we have used the fact that T b 2 = T s 2 = T b • T s = CV = 4/3 with Cf the usual Casimir 
of the fundamental representation of SU(3), E g is the soft-gluon energy and t is a unitary angular 
variable 

1 — cos-# 

t = (26) 



3 For simplicity's sake, we compute the square in Feynman (i.e. covariant) gauge, but the result is gauge 
invariant because the eikonal current is conserved: k fl J fJ '(k) = 0. 



5 



with d the emission angle of the gluon with respect to the strange quark (in the soft limit there 
is no recoil and the direction of the strange quark is not modified by gluon emission). The 
mass-correction parameter is given in this case by 

,2 



r = 



rn 



2 • 



(27) 



Note that in the massless limit for the strange quark, the last term on the r.h.s. of eq. (25) 
identically vanishes. The ordinary quark velocity and the Lorentz factor are given respectively 
by: 



v = 



El 



1 — r 
1 + r' 



1 = 



1 + r 



(28) 



By taking the ultra-relativistic limit r < 1, we obtain: 



Cf 



t + r 



1 



- 1 - 



{t + rf 



-g^J^k)r{k)* = 4na 

We now introduce the approximations 

1 _ 6{t-r) 
t + r ~ ' ~ ' (t + r) 2 

to extract the leading behavior of the matrix element squared for r <^ 1 4 . We then obtain 



(29) 



(30) 



C r 



9{t-r) . r0(t-r) 



t 



t 2 



C 

4na-^9(t-r) 



1 r 

i- l -7- 



■ (32) 



On the last member we have multiplied the term independent from t (related to an isotropic 
soft gluon emission off the beauty quark at rest) by the 6 function, as this only introduces 0{r) 
terms and allows a factorization of this "dynamical" constraint. 

The approximations (30) actually constitute the dead-cone approximation, which is the fol- 
lowing limitation on lower emission angles of the gluon with respect to the strange quark: 

m* 



Ei°y 



(33) 



where E^ is the strange quark energy in lowest order, which can be identified with the jet 
energy in the soft limit 

(34) 



E 



Xs 



E^ = ^(l + r), 



*Eqs. (30) are intended to hold in an integral way: 



dt 



1 dt 



+ O(r); 



dt- 



(t + rf 



f 1 1 

J dt p + 0(r). 



(31) 
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and m h is the beauty (pole) mass. We then obtain for the angular variable the corresponding 
limitation 



t ^> t m [ri 



1 — COS fir 



(35) 



An analogous simplification can be made for the kinematical constraint of given jet mass, i.e. of 
fixed (cfr eq. (11)) 

y - m \ s ~ m l (36) 



m b ~ m 1 



with m 2 Xs = (p s + p g ) 2 - Let us note that y is a unitary variable equal to zero in the Born 
kinematics point m Xs = m s- For a jet containing one soft gluon: 



y = (l-r)u 



t + 



1 — r 



= u(t + r) ~ cot 



for t > r, 



where 



2E„ 



m b (l - r) 



(37) 



(38) 



is the normalized gluon energy. 



To sum up: the soft-limit mass distribution for a jet initiated by a massive quark can be 
obtained from the corresponding massless formula by simply adding the term explicitly propor- 
tional to m 2 s and imposing the dead cone effect: 



1 dTW 



r dy 



soft 



f du [ 

JO Jr 



dt 



Aia D^a Dira 



The first-order coefficients read [5]: 



A l = 



cut 



C F 



+ 



t 2 



[5(y-ut) - 5(y)}. 



(39) 



D 1 = 



(40) 



7T 7T 

The first term proportional to D\ on the r.h.s. of eq. (39) is related to soft emission off the 
initial beauty quark, in its rest frame, while the second one is related to soft emission off the 
fast-moving strange quark. 

Eq. (39) misses the contribution from hard collinear gluon emission. We assume that the 
latter can be obtained from the massless one [5] by simply imposing the dead cone effect. The 
complete one-gluon distribution therefore reads: 



r dy 



did / dt 

Jr 



A\ol Dia D\ar B\ol 



LO t 



+ 



CO 



+ 



where 



B 1 = - 



cut 2 
3Cf 

4 7T 



+ 



t 



[5(y-ut) - 5(y)}, (41) 



(42) 
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A check of our assumption is provided by the comparison with the fixed order Feynman diagram 
computation of the spectrum given in sec. 2.2 5 . The event fraction is given by: 

rv i fir 

R{y) = I f ^dy' = l + J2«(y) + ..., (43) 

where: ^ 

R w (y) = - [ du f dt 

J0 Jr 

Note that the formula above has the correct end-point value: R(l) = 1. 



H 1 1 — 6(ut — y). 



lo t 



LO 



LO 



t 2 



(44) 



2.1.1 Leading Order 

Let us start by considering the leading-order term proportional to the coefficient A x in eq. (41). 
To include leading higher-order effects, we replace the tree-level coupling with the running cou- 
pling evaluated at the gluon transverse momentum squared [7] (see later): 

a -> a(k]_), (45) 

where the gluon transverse momentum squared is defined as [2]: 

k\ = [l-rfmluh ~ E 2 g d 2 ioid < 1. (46) 

By integrating the leading-order term on the r.h.s. of eq. (41) and expressing the result in terms 
of the gluon transverse momentum, we obtain: 

f X dLO f 1 dt A ,,,,,, . 1 f<V "Mh yM dk 2 

/ — / -A 1 a(kl)5(y-L0t) = - ± Al a(k 2 ± ) 

Jo LO Jr t V Jm*v* k l 

i r<y 2 /(y+^ dk 2 ± A /l2 , 

y Jmly> k{ 



(47) 



where on the last member a smooth interpolation between the two regions specified by the 
minimum has been considered according to the approximation 6 max[y,r] ~ y + r. Eq. (47) 
reduces in the frozen coupling case to — A x a / y\og{y + r) . As anticipated in the introduction, we 
have therefore two different regions: 



5 A more rigorous and direct derivation can be obtained by using the massive splitting functions (see for 
example [6]). 

6 In general, step approximations are more convenient to obtain simple analytical results while continuous 
functions are best suited for numerical purposes. 
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I. high jet mass 



- compared to quark mass m s : 



y > r. (48) 

Up to the logarithmic accuracy we are interested in, one can extend this region up to 
y > r. Eq. (47) reduces to the massless case: 

1 f m b y dk 2 

' ' .l ; n (// J ). (49) 

iv " 



In this region, the jet mass is so large that, at the logarithmic level, no effect is left of 
m s 7^ 0. To first order, eq. (49) becomes — Aia/y\ogy; 

2. low jet mass: 

y < r. (50) 

Analogously to the previous case, one can actually extend this region up to y < r. Eq. (47) 
specializes to: 

1 fm 2 b y 2 /r ji 2 

/„ ^.a(*i). (51) 

It is worth noting that the effect of a non zero mass, m s 7^ 0, is a restriction on the upper 
gluon transverse momenta. To first order, eq. (51) becomes — A\ajy logr. In this low jet- 
mass region, in agreement with physical intuition, the s quark mass m s completely screens 
the collinear singularity and produces a logarithm of the quark mass m s normalized to the 
relevant hard scale m^. 

2.1.2 Subleading effects 

Let us discuss in this section the derivation of the sub-leading effects in eq. (41): 

1. the first term proportional to Di, related to soft emission off the initial heavy parton at 
rest, reads: 

f — j dtD ia (kl)5(y-ujt) ~ D ia (m 2 b y 2 ) -. (52) 

J0 ^ Jr V 

This term is not modified with respect to the massless case, again as expected on the basis 
of physical intuition; 

2. the second term proportional to Di, related to soft emission off the massive and fast-moving 
strange quark, reads: 



bV » dA- 1 



y + r J \y y + r 



This term therefore vanishes in the high jet-mass region y ^> r, where we recover the 
massless case (in which this term was absent from the very beginning); 

3. the term proportional to B ± , related to hard collinear emission off the strange quark, reads: 



C du £ d j Bl a(kl) 6(y -ut)c Bl a (^) -L- c B ± a (m 2 y) °^Ljl, (54) 



where we have neglected small terms beyond the logarithmic accuracy. This term is regu- 
lated by the non-vanishing strange mass and, as expected, we re-obtain the massless case 
with the infrared cutoff r on y. 

By collecting the various contributions computed in the previous sections, we obtain for the 
0(a) distribution: 

1^ = - Aj^±A + 2D ia l - + ( Bl - D x )a—^—. (55) 
r dy y y y + r 

Within logarithmic accuracy, we can make a sharp approximation to obtain: 

1 dT (1) { -A 1 a\ogrl + 2D 1 a\, for y < r; 

r dy \ -A ia *fL + + for y > r. 

The main effects of virtual corrections are included, as usual, by replacing the above functions 
with corresponding plus distributions. The event fraction reads in the high-mass region: 

-— dy' = 1 - -A.alog'y + (B 1 + D 1 )a\ogy (y > r). (57) 

By imposing that the event fraction is continuous across y = r, 7 

R( y -> r - ) = -> r + ) = 1 - log 2 r + (Bi + £>i) a logr, (58) 

we obtain the logr terms in the low-mass region: 

R(y) = 1 — Ai a logy logr H -a log 2 r + 2D 1 ailogy + (-£>i — Di)o;logr (y < r). (59) 

2 



7 That is equivalent to using the plus distributions. 
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2.2 Check with fixed-order computation 



We compare in this section our eq. (43) with the 0(a) Feynman diagram computation of ref. [8], 
where the decay spectrum has been computed by retaining the non-vanishing strange-quark 
mass. The relevant contribution is the one coming from the magnetic penguin operator 

° 7 = T^ m b,MS^)sa^bF^. (60) 



m b,Ws(^ b ) * s ^ ne ^ mass in the MS scheme and fj,b = 0(m^) is the renormalization scale. By 
omitting non-logarithmic terms, the fixed-order (fo) distribution in eq. (30) of the first reference 
in [8] reads: 

y^L = . Aia ^il + 2Dia i + {Bl . Dl) i (61) 

r dy fo y y y + r 

The above formula exactly coincides with the first-order expression in eq. (55) derived in the 
previous section. This comparison confirms the validity of the inclusion of the coefficients of 
leading and subleading contributions as given in eq. (41). 



2.3 iV-space 



In order to resum the distribution to all orders in a, one has to transform it to moment iV-space. 
The invariant mass distribution in iV-space is defined as: 



N 



-f 

JO 



dy(l-y) 



N- 



r dy 



(62) 



The transform to N moment space indeed allows us to obtain an expression more suitable for the 
resummation of the logarithmic terms. As discussed in [9, 10] for the case of the jet mass event 
variable, the Mellin transform yields to a factorized expression of the phase space 8 . That gives 
rise to expressions for the jet mass distribution which are accurate at the next-to-next-to-leading 
order [9, 10]. The resulting expressions can be resummed in straightforward way. We will follow 
here the same procedure in order to resum the mass logarithms in eq. (41). The exponent in the 
resummed expression reproduces the usual structure, common to several resummed variables, 
as the one shown in refs. [2,9, 10]. In practice, one exponentiates the one-gluon distribution in 
N— space to account for multiple emission: 



N 



+ 



exp 



TV 



(63) 



8 The factorization of matrix elements follows instead from the dynamics of QCD. 
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Gluon branching, i.e. secondary emission, is taken into account by evaluating the running 
coupling at the gluon transverse momentum squared [7]: 



a — > a(k ± ), 

The "effective" one-gluon distribution therefore reads: 



(64) 



N 



[ dy(l- yf- 1 [ du [ dt 

Jo Jo Jr 



LOT UJ UJl z 



[8(y-ujt) - 5(y)}. (65) 



The virtual contributions can be "transferred" to the moment kernel as usual: 



TV 



f dy [{1 - y)"' 1 - 1] fdu f 

JO Jo Jr 



dt 



a(k]_) ha(fcl) ha(fc^ 



c<j t 



+ 



%-^). (66) 



In the following sections we perform the integrations of the various terms above. 
2.4 Resummation 

By putting all the pieces together, we obtain the jet-mass distribution for the heavy flavor decay 
in iV-space: 

T\ r /•! A,. { /•mjjj/min[l, 2 //r] ju2 

^ ~ exp/ ^[(l-^-l] / TT^i a + % - r)Si a (mgy) + 

1 JO !/ I JmW k ± 



+ D x a (m 2 b y 2 ) + 9(r - y)D 1 a (m 2 b y 2 /r) \. 



(67) 



A consistent next-to- leading logarithmic (NLL) resummation can be realized as in [2]. The first- 
order coefficients Ai, D\ and B 1 have been explicitly computed in the previous section. The 
remaining NLL terms are related to two-loop effects. In order to take them into account, in 
practice one has to include: 

1. the two-loop correction (oc in the QCD coupling a in the leading term A\ a; 

2. the two-loop correction to the first-order term A\ a by means of the replacement 

A x a -> A x a + A 2 a 2 . (68) 
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It is pretty well established that the latter function is universal - in the framework of 
effective theories, it is the well-known cusp anomalous dimension [11] 

r cusp (a) = r« p a + rf u ) p a 2 + • • • . (69) 



What about resummation in higher orders, i.e. NNLL and beyond? The coefficients A ± , A 2 , B x 
and D\ represent the lowest order terms of the functions: 9 

oo oo oo 

A( a ) = j2 A - a ^ B ( a ) = J2 Bna ^ D ^ = J2 DnaH - ( 7 °) 

n=l n=l n=l 

The functions A(a) and B(a) are related to small- angle emission only and therefore represent 
universal intra-jet properties, as confirmed by explicit higher-order computations [13]. The 
inclusion of the terms proportional to (3 2 , A 3 and B 2 is therefore rather safe — all these coefficients 
refer to universal ultraviolet or intra-jet quantities [14]. 

On the contrary, the function D(a), being related to soft emission at large angle, is in general 
a process-dependent inter-jet quantity. In the framework of fragmentation functions, the function 

DW(a) = D[ f) a + D ( 2 f) a 2 + ■■■ (71) 

has been originally introduced in [4], where it was called H(a); the first-order coefficient Hi 
was also explicitly computed — see the discussion in [15]. The D-function entering heavy flavor 
decays 

D«(a) = D[ h) a + D ( 2 h) a 2 + ■■■ (72) 
has however been shown to coincide with the former one to all orders [16]: 

D {f \a) = D w (a). (73) 

The function D^(a) with i — f,h refers basically to soft radiation emitted by a heavy parton 
with a small virtuality. In heavy flavor decays such as 

b^X s + ll (74) 

one considers a heavy flavor in its rest frame and looks at the final invariant mass distribution, 
while in heavy flavor fragmentation processes such as 

Z Q -> b + b , (75) 

9 A compilation of the known coefficients Ai, Bi and Di in our normalization, with references to the original 
papers, can be found in [12]. 
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one has a fast-moving beauty quark and looks at its energy. Since the function D(a) is the same 
in both cases, that seems to imply that the kinematical observable is irrelevant as long as soft- 
enhanced quantities are concerned. The only thing that matters is that of having a heavy quark 
close to its on-shell point point. The first order coefficients of the above functions D^\a) and 
D^ h \a) coincide with our D±. As discussed in the introduction, our D(a) is related to massive 
partons but it represents small-angle intra-jet corrections. On the basis of physical intuition, we 
therefore conjecture that this coincidence extends to higher orders, i.e. that our D(a) coincides 
with D^\a) or D^ h \a). From the above discussion it is clear that the inclusion of all the NNLL 
terms is not as rigorous as in the NLL case, the weakest point being the inclusion of the D 2 
term. An explicit check of our guess can be obtained by comparing our expanded resummation 
formula (given at the end of the paper) with an explicit (massive) two-loop computation, as 
soon as the latter becomes available. 

Beyond NNLL approximation, only the coefficient B 3 is known analitically [13]. We assume 
that the resummation formula keeps the same structure and that by including (unknown) higher- 
order terms to our expansions, 

A x a + A 2 a 2 + A 3 a 3 -> A(a), B x a + B 2 a 2 -> B(a), D x a + D 2 a 2 -> D(a), (76) 

we obtain the decay spectrum formally resummed to all orders: 

r„ f 1 A,, { rm 2 ymin[l,y/r] ii2 

- exp / ^ [(1 -„)"-!- 1] / ^A[a(kl)] + 9(y-r)B[a(m 2 y)] 

1 JO V [ Jm 2 b y 2 fc _L 

+ D [a (m 2 y 2 )} + 9(r - y) D [a (m 2 y 2 /r)} J. (77) 

In the limit r — > 0, we recover the well-known massless result. It is remarkable that the single 
logarithmic terms B and D have ^-functions of opposite arguments and therefore exclude each 
other: either a soft contribution is present, for a small jet-mass, or a collinear one is present, for 
a high jet mass. We now explicitly factor out the massless contribution in order to obtain the 
mass-correction factor 5n given in the introduction: 



r^v 

— = exp 

+ 6(r-y) 



f [d " yf- 1 ~ 1] ( C ^ [« (*i)] +B[a (mty)} + D [a (m 2 b y 2 )] 

Jo y [ Jm^y 2 K ±_ 



<V [a (k 2 ± )] -B[a (m 2 b y)] + D [a (m 2 b y 2 /r)] 



(78) 
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In the massive case, one has therefore the additional factor: 10 



S N {m 2 b ; ml) 



exp r£[(l-„)»-i-l]J - P" d -&A[c(kl)} + 

JO II I JmW/r fc ± 



B [a (m 2 b y)] + D [a (m 2 b y 2 /r)] 



(80) 



The integral above can be transformed back to unitary range by means of the rescaling v 
which gives: 



y/r, 



5 N (m 2 h ] m 



1 dv 



exp / — [(1 — rv 
'o v 



,N-1 



-i] - 



dk 2 



m 2 v 2 ft j 



A [a (hi)] + 



- B [a (m 2 s v)] + D[a (m 2 s v 2 )] \. (81) 



It is remarkable that, as a consequence of the rescaling, the hard scale Q (appearing in the 
limits of the transverse momentum integral as well as in the argument of the coupling in the 
single-logarithmic terms) changes from mj to m s . Let us also note that, in the last equation, 
the hard scale m b only enters through the variable r inside the moment kernel. 



To express the correction factor as a standard Mellin transform, we use the relation: 



(l-ryf- 1 - 1 = (l-yr N -V - 1 + O 



(82) 



Eq. (82) can be easily shown to be valid at the next-to-next-to-leading log N level, by using the 
relation [17] 



(i _ y) »-i_ i ^ _ * I „_i) + £2> 



n 



n 



7?/ 



n 



+ O 



N 3 LO, — 



where n = Ne 111 , 7^ = lim 



k— >oo 



0.577216' 



z(a) = Yl^=i^/ na i s ^ ne R-iemann Zeta function with z(2) = 7r 2 /6 = 1.64493 The final 
expression for the correction factor in iV-space therefore reads: 



(83) 

is the Euler constant, 



S N {m 2 b ; m 2 ) 



exp 



1 dv 



2, ,2 k 1 



■A [a (k\)] + 



10 A compact derivation of the correction factor 5n is obtained by splitting the angular integral in all the terms 
on the r.h.s. of eq. (66) but the term proportional to D\r, as: 



r-l i*l r-r 

I dt= dt- dt. 

Jr JO Jo 



(79) 



The first integral on the r.h.s. represents the massless case while the second one the mass corrections. The 
unitary range in the latter is restored by setting if = t/r. 
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- B [a (m 2 s v)] + D [a (m 2 s v 2 )] \. (84) 



3 General case 

This is the central section of the paper and contains general results about threshold resummation 
in processes with jets initiated by partons with a small mass compared to the hard scale. 

3.1 Mass Effects in a Jet 

We generalize in this section the resummation formula obtained in the previous section for the 
radiative b decay by simply noting that any reference to the particular process disappears in the 
correction factor Sn in eq. (84). Therefore we simply replace the beauty mass with the hard 
scale Q of the general process and m s with the mass m of the quark triggering the jet under 
consideration: 

rl r(JV-l)_i f r m 2 (l-z) J k 2 

5 N (Q 2 ;m 2 ) = exp / dz -I- [a (k 2 ± )} - B [a (m 2 (l - z))] + 

J0 Y - Z y Jm 2 (l-z) 2 fc ± 

+ D[a(m 2 (l-z) 2 )}^, (85) 

where r is now defined in eq. (13). A check of our generalization is provided by the comparison 
with the full 0(a) computation of the DIS cross section with a massive quark in the final state 
(see next section). 

As anticipated in the introduction, the jet mass distribution is an infrared safe quantity, 
i.e. it has a finite limit for vanishing quark mass, m — > 0. Our results are in agreement with 
this general fact in the following way: for r ^ 0, the distribution contains mass logarithms 
~ log r which are not power suppressed, but have support in the region of power-suppressed size 
y < r 11 . The 0(a) term in the correction factor in physical space 5(y; Q 2 ,m 2 ) can be obtained 
by subtracting out line 2 (the massless distribution) from line 1 (the massive distribution) in 

n Lct us note that logarithmic mass effects only occur for rN » 1 because, for rN <C 1, one can expand the 
exponent in (85) as 

z r ( N -V - 1 = r(iV-l)logz + \r 2 (N- l) 2 log 2 z + ••• , (86) 

obtaining power-suppressed effects of the usual form. The use of the resummation formula in the latter case 
however is not legitimate. 
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eq. (56): 

S(y; Q 2 ,m 2 ) = 5{y) - A x a 



9{r - y) 



\ogrJy_ 

y 



+ (Di - B x )a 



Ojr - y) 

y 



+ 0{a 2 ). (87) 



+ 



Let us note that 5 actually vanish for y > r (where r < 1). 



3.2 Check with DIS with massive final quark 

The first order corrections to the inclusive cross section have been computed in [18] for the charm 
production in charged- current DIS, 



(88) 



where (g) is a real or a virtual gluon. In this computation, the non zero charm mass has been 
retained while the (much smaller) strange quark mass has been neglected: 



m r 



m ^ 0; m s = 0. (89) 
Omitting non-logarithmic terms, the cross section given in eq. (40) of [18] can be written as: 



+ 



CfOL 



TV 



2 log(l — x) — log(l — \x) 



l-x 



Cpa 



+ 



1 C F a 1 



+ 



7i [1 - x] 4 7i [1-Xx] 



where 



A = ^ 



(90) 



(91) 



Q 2 + m 2 ' 

x = xb = Q 2 / (2p • q) is the Bjorken variable, Q 2 = —q 2 is the hard scale squared, q = — p V)1 
is the P^-boson momentum, /ip is the unit of mass of dimensional regularization (to become the 
factorization scale after pole subtraction), Pqq\x) is the leading-order (massless) q — * q splitting 
function in 4 dimensions, 



C f 



1 — x 



(92) 



and 



-z = - ~ 1e + log(47r), 



(93) 



with d = 4 — 2e the space-time dimension. For simplicity's sake, let us set from now on jip = 
Q. The infrared pole, of collinear nature, is absorbed into the quark non-singlet distribution 
function, which reads, in the MS scheme: 



F q/q (x; Q 2 ) = 5(1 -x) - 



a 1 
2tt i 



P^{x)+0{a 2 ). 



(94) 
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The next step is to factorize from the cross section in eq. (90) the (massless) coefficient function: 



log(l - x) 



1-x 



C DIS {x; Q 2 ) = 5(1 -x) + A l a 

The latter is obtained as the following convolution [2-4] 

C D is(x; Q 2 ) = A(x; Q 2 ) ® J(x; Q 2 ), 

where 



[1-x] 



+ O (a 2 ) . (95) 



+ 



A(x; Q 2 ) = 5(1 -x) + 2A ± a 



log(l — x) 
1-x 



+ O (a 2 ) , 



(96) 
(97) 



+ 



is the radiative factor related to the observed initial-state jet, produced by the massless s quark, 
while 



J(x; Q 2 ) = 5(1 -x) - A 1 a 



log(l — x) 
1-x 



+ O {a 2 ) 



(98) 



is the jet factor related to the unobserved final-state jet, initiated by the c quark, in the massless 
approximation. The cross section can therefore be written as: 



where: 



~ F(x; Q 2 ) <g) A(x; Q 2 ) ® J(x; Q 2 ) ® h(x; Q 2 , m 2 ), 

a ax 



h(x;Q 2 ,m 2 ) = 5(l-x)+A l a 



(99) 



! 


log(l — x) 




log(l — x + r) 




1 — X 


+ 


1 — X 



+ 



+ (D 1 - Bx) a 



[l-x] + [1-x + r] 

The "®" denotes a convolution and we have used the fact that 

1 



A 



1+r 



~ 1 — r 



and 



1 - A x = 1-x + r. 



(100) 

(101) 
(102) 



The last term h(x; Q 2 , m 2 ) on the r.h.s. of eq. (99) is in agreement, within logarithmic accuracy, 
with 5(y; Q 2 ,m 2 ) expanded to 0(a), in eq. (87), after setting x = 1 — y: 



h(x; Q 2 , m 2 ) ~ 5(1 - x; Q 2 , m 2 ). 



(103) 



The resummation of mass effects in the charged- current DIS cross section, for several values of 
the ratio m c /Q, has been investigated in detail in [19]. 
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3.3 Tower Expansion 



The universal mass-correction factor has the generalized exponential structure [2] 

5 N (Q 2 ; m 2 ) = eM« 2 ;™ 2 ), (104) 
where the exponent has a double expansion of the form: 

oo n+1 

F N {Q 2 ;m 2 ) = ^^F nfe a"log fc (iVr), (105) 

n=l k=l 

with F nk numerical coefficients. The exponent can be expanded in towers of logarithms as: 



F N (Q 2 ;m 2 ) = pd x {p) + X> nrf n +2 (p) 



ra=0 



pd 1 (p) + d 2 {p) + ad 3 (p) + a 2 d A (p) + 



where 



p = p a(p 2 ) log (Nr) 
and p = 0{m) is a renormalization scale of the order of the quark mass m. 



(106) 
(107) 



By truncating the above series expansion, one obtains a fixed-logarithmic approximation to the 
form factor 5^. The functions di(p), which represent the mass effects, can be obtained from the 
standard ones (?i(A) [12] by means of the replacements: 12 



A(a) 



A(a); B{a) 



2 2 

B(a); D(a) - D(a); log-^-^log^; A - p. (108) 



Q 2 



It is worth observing that mass effects induce a similar structure to the massless one, involving 
changes of sign of the collinear functions A and B, with the rescaling Q — > m. We then obtain: 



di(p) = 
d 2 (p) = 



A 1 

2/V 



1-2 p) log(l-2p)-2(l-p) log(l-p) 

A 2 



+ 



2Po 
AiPi 

4/?o 3 

A\lE_ 



El 



log(l-2p)-^log(l-p)- 



(109) 

log(:l - 2p) -2 1og(] p)\ + (110) 

+ 



2 log(l - 2 p) + log 2 (l - 2 p) - 4 log(l - p) - 2 log 2 (l - p) 

A l 



log(l-2p)-log(l-p) 



2/9 



log(l-2p)-2 1og(l-p) log 



12 All these functions contain in principle the over-all factor 9 (N — 1/r), coming from the step approximation 
of the moment kernel, which avoids modifications for small N of the massless behaviour, in agreement with 
physical intuition. Analytic continuation to the complex TV-plane is made by omitting such factor. 
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For the NNLO function we obtain: 



+ 



+ 
+ 



+ 
+ 



+ 



+ 



Do 



A) l-2p 



2p , Z?i /3i 



+ 



#2 P R P 

A) 1 - P 1 - P 



l-2p 2 A/ 
+ 



2p t log(l-2p) 



+ 



l-2p l-2p 



Bilk 
A) 2 



P , lo g(l - p) 
1-P 



1-P 



+ 



+ 



A, 



2 A/ 



P 



l-2p 1-p 



+ 



A 2 IE 


2p P 


A) 


1 


-2p 1-pJ 




3 p 3 p 


2 A) 3 [ 


1 - 


-2p 1-p 


A^e 2 




4p P 


2 


1 


-2p 1-p 


a^ 2 r 


2p 2p 


4 A) 3 [ 


1 - 


-2p 1-p 


Ai A 7i? 


r 2p p 


A) 2 




l-2p 1- 


AiA 2 




P P 



+ 



+ 



1 -2p 

M 7T 2 

12 



4p 



1-p 
P 



l-2p 1-p 



2 /3o 4 



l-2p 1-p 



log(l-2p) _ log(l-p) 
p l-2p 1-p 

-,og(l-2 rt + 1 ^ 1 - 2 ^ 



+ 
+ 

+ 
+ 



log 2 (l-2p) 



2 (l-2p) 

m 



+ 2 1og(l-p) 
1 B 1 p 



l-2p 
2 1og(l-p) log 2 (l-p) 



A) i-2p 



A 

A) 3 



2p 



P 



p 

Pol- p 

2 



1-p 

2 A, 



1-p 
P 



+ 



P 



l-2p 1-p 



+ 



i P 

log^ + 



l-2p 
P 



1-p 
P 



i P 

log At + 



- log(l-p)- 



2p 1-P 
P log(l - p) 



| log(l-2p) | plog(l-2p) | 



2p 



2 A) 



2p 2 



2p 



1-p 

P 2 
1-P. 



, P 



i 2 P 

log 



(HI) 



The coefficients A °f the QCD /3-function in our normalization have been given in [12]. 
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3.4 Inverse Mellin Transform 



The mass- correct ion factor in physical space is obtained by means of an inverse Mellin transform 
of 5n- 

« (»;«»,».») = (!-») | { J T°^(i-y)-^(0", m ')}. (H2) 

where c is a (real) constant chosen in such a way that the integration contour lies to the right 
of all the singularities of S N . By defining 

5 Nr = 5 N (113) 
and changing variable from N to u = Nr, we obtain: 

$ (</; q 2 ,"! 2 ) = £ {£j^ ^-yf'TUQ^)}, (114) 

where we have omitted the power correction y multiplying the derivative with respect to y. The 
correction factor in physical space is therefore the inverse Mellin transform of S u with respect to 
(1 - yf/ r 13 . We can therefore use the results in [17] to obtain the correction factor in physical 
space in NNLL approximation: 

5 (y; Q\ m 2 ) = A {%) A (y; Q 2 , m 2 ) } (115) 



dy 

where: 



e Ld 1 (r) + d 2 (r) 



A ^ Q2 ^ = w^\ SA - (116) 

We have defined 

Mr) = rKWl = d^ + rd^r). (117) 
ar 

8 A is a NNLL correction factor which can be set equal to one in NLL: 

5&NLL = I- (118) 

Its NNLL expression reads: 



5A = (119) 



with 



S = e ad ^ll +/3oa4(T)V[l-/ii(r)] + ^(3 a h[(r) {^j 2 [1 - d^r)] - ^j' [1 - d l (r)] 



(120) 



13 Let us note that, in the large N limit we are interested in (see eq. (82)) rN 3> 1, we can make the approxi- 
mation r(N — 1) ~ rN — 1. 
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T(x) is the Euler Gamma function and 

^j(x) = ^ogT(x) (121) 

is the digamma function. Finally: 

t = (3 aL (122) 

and 

L = -logfl-a-!/) 1 ^. (123) 
It is convenient to approximate the argument of the logarithm by an expansion for y <^ r: 

(i-y) 1/r = i- Z + o(£). (124) 



Note that the r.h.s. is positive only for y < r, implying that the linearization above shrinks the 
domain of y from the unitary interval (0, 1) to the much smaller interval (0,r). 

To summarize, we have the final result: 

5 (y; Q 2 , m 2 ) = L { d (y) A (y; Q 2 , m 2 ) } (125) 
where A (y; Q 2 , m 2 ) is given in eq. (116) and 



r 



L = 9{r-y) log-. (126) 

y 

In agreement with the observation above, we have limited the domain to y < r with a ^-function. 



3.4.1 Coefficients of the mass logarithms 

The mass correction factor in physical space A(y;Q 2 ,m 2 ) also has a generalized exponential 
structure: 

A(y;Q 2 ,m 2 ) = e H ^ Q "^) , (127) 
where the exponent has a double expansion of the form: 

oo n+1 

H (y; Q 2 , m 2 ) =6(r-y)J2Yl H ^ ^ lo S* ( 128 ) 

n=i k=i y 

with H n k numerical coefficients. By expanding the r.h.s. of eq. (116) up to third order, one 
obtains the following expressions for the logarithmic coefficients: 

#12 = \A V) (129) 
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#11 

#22 
#21 
#34 

#33 

#32 



B 1 - D 1 ; 
1 



lA 2 + ±p (B 1 -2D 1 )-±Alz(2); 
B 2 -D 2 - A, (B 1 - D x ) z(2) - A 2 lZ (3); 

T2 A ^ 

A 2 (3 + \a x ^ + (B 1 -AD,)- ^(2) " ^(3); 
l -A z + (3 (B 2 -2D 2 ) + ^(B l -2 D 1 ) - A, A 2 z{2) - (5^-7 D,) z{2) + 

Al %(4)-Mi!|^)_V(Bi-^iM3), 



(130) 
(131) 

(132) 
(133) 
(134) 

(135) 



1 ' 2 ' - ™ < 136 > 

where ,(3) = 1.20206 • • • and z(A) = ?r 4 /90 = 1.08232 • • • . Note that the leading coefficients H 23 
and H34 involve products of the one-loop coefficients A\ and (3q only. The explicit expressions 
of the coefficients read: 



#12 
#11 

#23 
#22 
#21 



2tt' 
4tt' 

C f fllC a 



Air 2 



6 



3 



C 2 F z(2) C F 
2ir 2 4tt 2 



Ca[^-z(2) 



C f 



#34 
#33 



4tt 2 
48tt 3 



5,(3) ) 
239 2,(2) \ 



108 



+ C A 



35rif 
36 



-3595 5,(2 _ ,A 

+ — ^ + 19, 3 

216 3 V 7 



(137) 
(138) 
(139) 

(140) 



C{ 8A7C A 77n f } | 7n? 



36 



9 



9 



(141) 
(142) 



<%*(3) Cj. 
3tt 3 4tt 3 



11CW2) , 1 

^+n/(-4+«(2)) 



C f 



An 3 



lln2 f +c \( ini m(2) ^ 



36 



V 144 



6 



) 



+ C A n } -4 + 



,(2) 



(143) 



#32 = 



( * [z(A) + ,(3)] + 



47T 3 
471-3 



-67rif 61nfz(2) 



48 



36 



+ 5n f z(3) + C A (-g - + 2 ,(2) 2 - 22,(3)) 



+ 



47T 3 



-3n 2 n 2 ,(2) ^ -8855 145,(2) 11,(4) 319,(3) ^ 

~^ + ^ _ + [~Wa 3^ + ^A~ + ^2~) + 
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(144) 



where Ca = N c = 3 and rif is the number of active quark flavors. 

4 Conclusions 

We have resummed the invariant mass distribution of hadronic jets initiated by massive quarks 
in next-to-leading logarithmic approximation. The resummation has later been extended in 
heuristic way to the next-to-next-to-leading logarithmic approximation. The mass effects have 
been relegated into a universal factor 5n which takes into account the well-known dead-cone 
effect and soft radiation characteristic of massive partons. 5n contains the same resummation 
functions which are encountered in standard threshold resummation. Mass corrections produce a 
universal intra-jet structure in agreement with one's physical intuition: only small angle partons 
emissions are involved, which can be ascribed to a specific jet in the event. It is remarkable that 
the coefficients of the mass logarithms are simply connected to those found in massless processes, 
at low order as well as in higher order of perturbation theory. A similar situation is found in the 
fragmentation of heavy quarks [4] . Our formulae have been checked against explicit first-order 
computations: the radiative b decay b — > S7 with m s 7^ and DIS + s — > c + /i with m c 7^ 0, 
finding complete agreement. 

Mass effects, i.e. effects related to non vanishing parton masses, often play a significant role 
in jet physics at the quantitative level [20]. It may be worth to cite just a few applications of 
our results. 

Semi-inclusive £>-decays 

B -»• X c + I + v { (145) 

are largely affected by the non- vanishing charm quark mass as r >rn 2 c /ml ~ 0.1. Semileptonic 
b — > c decays allow the extraction of the CKM matrix element V c b and constitute the main 
background to the b — > ulv decays, which are used for the extraction of the CKM matrix 
element V u b- The inclusion of charm mass effects is needed to have a better understanding, for 
example, of the charged lepton spectrum or the invariant hadron mass distribution, which have 
recently been measured with great accuracy in [21-26]. The resummed formula which we have 
obtained can be combined with the full 0(as) triple differential distribution for (145) recently 
obtained in [27, 28]. An additional complication in this case stems from the fact that the charm 
quark velocity, or equivalently the parameter r, is not fixed in the tree-level process b — > civ. 
That is because the hard scale, set by the total hadron energy in the final state, is not fixed. 
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Since the exclusive channels B — > Dlv and B — > D*lv constitute a large fraction of the total rate 
(145), an 0(50%) [29], the use of perturbation theory can be questioned. We believe however 
that a perturbative computation can provide quantitative informations on the decay (145) to be 
compared with other models. 

An accurate computation of shape variable distributions in e + e~ annihilations at the Z° 
pole and below — such as thrust, heavy jet mass, C-parameter, etc. - - asks for the inclusion 
of the beauty mass effects very close to the two-jet region. The mass correction parameter 
r w 4ml /s ~ 0.1 for y/s = 30 GeV, while it is smaller by an order of magnitude at the Z° pole. 
At a future e + e~ linear collider of center-of-mass energy of 500 GeV, mass effects in top pair 
production will be controlled by r ps 0.5. 

The invariant hadron mass distribution in semileptonic top decays, 

t -> X b + W (146) 
is affected by the non-vanishing beauty mass close to threshold, i.e. for mxb > mb- 

In ref. [30] a disagreement has been found in the (massless) evolution of charm fragmentation 
data from 10 to 91 GeV. We argue that the inclusion of charm mass effects, which should be 
significant for N > (m c /5)~ 2 w 10, could improve the accuracy of the perturbative computations 
and eventually solve this problem. 

To sum up, mass effects in a jet can be included, within logarithmic accuracy, by the universal 
factor in eq. (18) multiplying the massless jet function. 
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